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J> ■ Abstract 

' All complex Hadamard matrices in dimensions two to five are known. We use this 

C^) . fact to derive all inequivalent sets of mutually unbiased (MU) bases in low dimensions. 

We find a three-parameter family of triples of MU bases in dimension four and two 
f — , . inequivalent classes of MU triples in dimension five. We confirm that the complete sets 

of (d + 1) MU bases are unique (up to equivalence) in dimensions below six, using only 
C^S . elementary arguments for d less than five. 

> ■ 

;rj ■ 1 Introduction 

rN ' 

Position and momentum of a classical non-relativistic particle are intimately linked since 



the momentum variable generates spatial translations. Mathematically, this important 
relation is embodied in the structure of the Galilei group. It turns out to be even more 
fundamental for a quantum mechanical particle, where it takes the form of the commutation 
relation of its position and momentum operators. 

The associated Heisenberg-Weyl group of phase-space translations continues to be rel- 
evant for quantum systems with only a finite number of orthogonal states, providing a 
basis of the space C d . For each dimension d > 2, there is a set of unitary operators which 
give rise to a discrete equivalent of the Heisenberg-Weyl group [I] . Physicists would expect 
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the state spaces C to be structurally identical, at least with respect to properties closely 
related to the Heisenberg-Weyl group. It thus comes as a surprise that the Heisenberg- 
Weyl group allows one to construct (d + 1) so-called mutually unbiased (MU) bases of the 
space C d if d is the power of a prime number [21 [3] while the construction fails in all other 
'composite' dimensions. No other successful method to construct (d + 1) MU bases in all 
dimensions is known [U (5J [6] . 

Given (d + 1) orthonormal bases in the space C d , they are mutually unbiased if the 
moduli of the scalar products among the d(d + 1) basis vectors take these values: 



where b, V = 0, 1, . . . , d. Such complete sets of MU bases are ideally suited to reconstruct 
quantum states [3] while sets of up to {d + 1) MU bases have applications in quantum 
cryptography [8] and in the solution of the Mean King's problem [9], for example. 

The methods to construct complete sets of MU bases typically deal with all prime 
or prime-power dimensions simultaneously. They either make use of the Heisenberg-Weyl 
group [10] , exploit identities from number theory j3[ [TTJ , or they are couched in the language 
of finite fields [3j [12] . All these methods are constructive and effectively lead to the same 
bases. The existence of other, inequivalent sets of complete MU bases remains unclear. 

In this paper, we choose a different method to study MU bases in dimensions two to 
five. It allows us to directly conclude that the corresponding complete sets of (d + 1) MU 
bases are unique (up to some irrelevant equivalence, cf. below). What is more, we are also 
able to exhaustively list all inequivalent classes of d or less MU bases. This approach is 
attractive because it uses elementary methods only, except in dimension five. In order to 
present the details of this method we need to briefly discuss the relation between MU bases 
and complex Hadamard matrices. 

MU bases and complex Hadamard matrices 

Each MU basis in the space C d consists of d orthogonal unit vectors which, collectively, 
will be thought of as a unitary d x d matrix. Two (or more) MU bases thus correspond to 
two (or more) unitary matrices, one of which can always be mapped to the identity / of 
the space C d , using an overall unitary transformation. It then follows from the conditions 
([T]) that the remaining unitary matrices must be complex Hadamard matrices: the moduli 
of all their matrix elements equal \j\fd. This representation of MU bases links their 
classification to the classification of complex Hadamard matrices |13t [Tl] . 

It is, for example, possible to list all pairs of MU bases {/, H} in C d once all complex 
Hadamard matrices are known. Using this observation as a starting point, we will extend 
the classification from pairs to sets of r < (d + 1) MU bases. As complex Hadamard 
matrices have been classified for d < 5, we expect to obtain an exhaustive list of sets of r 
MU bases in these low dimensions. 




if b = b' 
ifb^V 



(1) 
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The approach we take is inspired by recent works aimed at dimension six [151 [T6l [T7] . 
Starting with a pair of MU bases {/, H} we will find all vectors v which are MU to / and 
H. Only these vectors represent candidates to form additional bases and by analysing their 
inner products we will be able to obtain all MU bases in low dimensions. 

In dimension six it is not trivial to determine all candidate vectors. It becomes necessary 
to use Grobner bases |16j or to discretise the underlying space |17j . For d < 5, however, 
we find that elementary properties of the complex plane are sufficient to solve (most of) 
the relevant equations in closed form. 

The task to find all MU bases is complicated by the fact that, actually, many sets of 
apparently different MU bases are identical to each other. For the desired classification, it 
is sufficient to enumerate all dephased sets of (r + 1) MU bases. This standard form |13| is 
given by 

{I,H 1: ...,H r }, re{l,...,4, (2) 

where I is the identity in C d and the other matrices are complex Hadamard matrices of a 
particular form: the components of the first column of the matrix Hi are given by 1/Vd, 
and the first row of each Hadamard matrix H p , 1 < p < r has entries \j\fd only (see Eqs. 
(133p and (|35h for an explicit example with d = 5). The possibility of dephasing is based on 
the notion of equivalence classes for MU bases, explained in more detail in Appendix A. 
Roughly speaking, two MU bases are equivalent to each other if one can be obtained from 
the other by changing the overall phases of individual vectors and by permuting them. 

The results of this paper have been arranged as follows. In Sec. [2] we deal with 
dimensions two and three. The complete list of sets of MU bases in dimension four is 
derived in Sec. [3j Then, all sets of MU bases of C 5 are constructed, and in Sec. [5] we 
summarize and discuss our results. 



2 Dimensions d = 2 and d = 3 

In this section, we construct all sets of MU bases in dimensions two and three using only 
simple properties of the complex plane. The direct approach to construct all MU bases for 
d = 4 in Sec. [3] will be based on similar arguments. 

2.1 Dimension d = 2 

The matrices consisting of the eigenvectors of the Heisenberg-Weyl operators form a set 
of three MU bases in dimension two which are unique up to the equivalences specified in 
Appendix A. We present a simple proof of this well-known fact. 

Let us begin by noting that there is only one dephased complex Hadamard matrix in 
d = 2 (up to equivalences), the discrete (2 x 2) Fourier matrix 
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A vector v G C 2 is MU to the standard basis I (constructed from the eigenstates of the 
z-component of a spin 1/2) if its components have modulus 1/y/d. Applying the transfor- 
mation given in Eq. (|54p . the dephased form of such a a vector reads v = (1, e ia ) T / y2, 
with a real parameter a G [0, 27r]. The vector v is MU to the columns of F 2 if the phase a 
satisfies two conditions, 

\l±e ia \=V2. (4) 

These equations hold simultaneously only if e ta = dbi. Thus, there are only two vectors 
which are MU to both I and F 2 , given by v± = (\,±i) T Jy/2. Since this is a pair of 
orthogonal vectors, they form a Hadamard matrix H 2 = and, therefore, the three 

sets 

{I},{I,F 2 },{I,F 2 ,H 2 } (5) 
represent all (equivalence classes of) one, two, or three MU bases in dimension two. 



2.2 Dimension d 



In dimension three there is also only one dephased complex Hadamard matrix up to equiv- 
alence. It is given by the (3 x 3) discrete Fourier matrix 



(6) 



defining uj = e 27 ™/ 3 . Again, we search for dephased vectors v = (1, e ia , e l ^) T /\/3, < 
a, (3 < 2ir , which are MU with respect to the matrix F$. This leads to the following three 
conditions 




ifi 



l + e ia + e 



1 + ue ia + uj 2 e i0 



1 + u 2 e ia + ue i0 



V3 : 

Vs. 
Vs. 



(7) 



Removing an overall factor of e ia / 2 , they can be rewritten 



. a 
C + cos - 



. a , 2ir\ 
C + cos(- ±T j 



^3 
2 

Vs 
2 



(8) 



where = a / 2 ). By considering a plot in the complex plane, Fig. [H we see that 
these three equations hold simultaneously only if two of the cosine terms are equal. This 
implies that the only possible values of the parameter a are 0, 7r/3, or 2n/3, leading to the 
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requirement ±1/2 = cos/3. Consequently, the Eqs. ([7]) have exactly six solutions which 
give rise to vectors 



Vl oc 



V4 OC 




Examining their inner products shows that there is only one way to arrange them (after 
normalization) into two orthonormal bases, namely H$ = (vi\v2\v 3 ) and H$ = (v^vslve). 
It is useful to note that one can write 

ijW = DF 3 and = D 2 F 3 , (10) 

where D = diag(l, u, uj) is a diagonal unitary matrix with entries identical to the com- 
ponents of the vector v\, i.e. the first column of JEZ3 . The triples obtained from adding 
either or to the pair {/, F 3 } are equivalent, 

{I,F 3 ,HP} ~ {D 2 ID,D 2 F 3 ,D 2 H { 3 1) } = {I,h£\f 3 } ~ {I,F 3 ,H^}, (11) 

as follows from first applying the unitary D 2 globally from the left, rephasing the first basis 
with D~ 2 = D, and finally rearranging the last two bases. We therefore conclude that the 
sets constitute a complete classification of all sets of MU bases in dimension d = 3. 



3 Dimension d = 4 



In dimension d = 4, a one-parameter family of complex Hadamard matrices exists, 
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x e [o, 7r] 



(12) 



It is elementary to show that all 4 x 4 complex Hadamard matrices are equivalent to a 
member of the family F±{x). When x = 0, the resulting matrix is the discrete Fourier 
transform F4 on the space C 4 , with matrix elements given by oji k , j, k = . . . 3, oj = i. The 
matrix F^(tt/2) is equivalent to a direct product of the matrix Fi with itself while for other 
values of x it can be written as a Hadamard product of F4 with an x-dependent matrix. 
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Im z 




Figure 1: Plot of Eqs. ([8]) in the complex z-plane. The real numbers c\ and C2 each 
represent one of the three numbers cos (a/2) and cos (a/2 ± 27r/3); it follows that at least 
two of these three expressions must be equal. 

3.1 Constructing vectors MU to F±{x) 

After dephasing, any vector MU to the standard basis takes the form v = (l,e ia , e^ 3 , 
e^') T /2 where < a',/3',7' < 2vr. For convenience, we will use an enphased variant of 
v. Multiplying through by the phase factor e~ ia I 2 and defining a = a'/2 S [0, 7r], f3 = 
ft 1 — a'/2 £ [0, 2tt], and similarly for 7, we consider the parametrization v = (e~ ia , e ia , e 1 ^ , 
e* 7 ) T /2 instead. The conditions for u(a,/3,7) to be MU to the columns of -^4(2:) lead to 
four equations, 

|cosa±£+| = 1, (13) 
|sina±e-"'C_| = 1, (14) 

where complex numbers C± = ( e *^ =t e* 7 )/2 have been introduced. We will now construct 
all solutions of these equations as a function of the value of x. We treat the cases (i) a = 0, 
(ii) a = 7r/2, and (hi) a/0,a/ 7r/2 separately since the Eqs. (|13p and fll4f) take different 
forms for these values. 

(i): a = 

Eqs. (I14p simplify to the pair |±e _llE C-| = 1 > which only hold simultaneously if |C-| = 1 
or e n = —e l/3 , implying that £ + = so that Eqs. (fT3|) are satisfied automatically. Thus, 
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solutions exist for any value of x whenever /3 = j + ir mod2-7r, and the resulting vectors can 
be written as v(j3) = fl, 1, e*^, — e'*) /2, with j3 G [0,2-7r]. It will be convenient to divide 
this family of states into two sets, 



/ 1 \ 



hi(y) 



\ 



o'.V 



( 1 \ 



, h 2 (y') 



< y,y' < it, 



(15) 



introducing y = (3 and y 1 = it + (3. 
(ii): a = vr/2 

Eqs. f)13|) and (|14p now reverse their roles: the conditions | ± £ + | = 1 require e* 7 = e*^, 
with (I13p being satisfied since |£_| = 1 follows immediately. Hence, there is another one- 
parameter family of mutually unbiased vectors for all values of x if f3 = 7. This family 
can be written as v(ip) = (l,—l,e itf ,e l(p ) T /2, tp £ [0, 2ir], after dephasing and absorbing a 
factor of i in the definition of the phase, ip = tt/2 + /3. Again, we express these solutions 
as a set of pairs, 



h 3 (z) = - 



( 1 \ 
-1 

e iz 



( 1 \ 



h A {z') 



< Z, z' < 7T, 



(16) 



where z = <p and 2/ = 7T + </?. 
(iii) a ^ 0,a ^ vr/2 

A plot in the complex plane (see Fig. [2]) reveals that one must have = ±i sin a if Eqs. 
(fT3|) are to hold with cos a^O. Thus, the real part of C+ vanishes, 



cos /3 + cos 7 = 



(17) 



with 7 = 7T — /3 mod 2ir being the only acceptable solution: the other solution, 7 = tt + 
(3 mod27r leads to = C+ = ±isina, producing a contradiction since a ^ 0. Thus, using 
7 = 7r — /3 mod27r, we obtain £ + = isin/3 find the following relation between a and /3: 



± sin a = sin /3 . 



(18) 



Similarly, Eqs. ()14[) for sin imply that e lx Q- = ±i cos a. Using 7 = 7r — /3 mod 27r 

in the definition of we find = cos /3, so that 



i(± cos q + sin x cos /3) = cos x cos /3. 



(19) 
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The right-hand-side of this equation only vanishes if x = tt/2: both f3 = tt/2 and (3 = 3ir/2 
would, according to (|18p . require a = tt/2 which we currently exclude. Therefore, solutions 
to Eqs. ()13ll4p with a / or a / tt/2 only exist for x = tt/2 if a second relation between 
a and /3 holds, 



± COS Q = COS P . 



(20) 



The form of the additional MU vectors is determined by Eqs. f)18|) and (|20j) which have 
four solutions. First, for /3 = a we obtain MU vectors of the form (e~ ta , e ia , e ia , —e m ) T /2 
or (l, e 2lQ , e 2lQ , — l) T /2 after dephasing. Splitting this family into two subsets as before, 
we find 



I 1 \ 



2 



e 

V -i J 



( 1 \ 



,k 2 =l 

2 



< t, t' < IT . 



(21) 



V -i J 

Similarly, the choice /3 = it + a mod27r leads to two sets of dephased MU vectors, 



1 



/ 1 \ 



V -i J 



( 1 \ 



n iv! 



, < U, U < 7T . 



(22) 



V -i / 



Next, when proceeding in an entirely analogous manner for the remaining two choices 
/3 = 7r — a mod27r and f3 = 2tx — a mod27r, we obtain the following four families of 
dephased vectors MU to F&(k/2), 



31 
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, J3=2 




1 




' H= 2 




1 
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e ir ) 




V -e ir ' 
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_ e is 
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K 


e is> 


) 



(23) 



with < r, r', s, s' < n. 



3.2 Forming MU bases 

Knowing all vectors that are MU to both the identity and F4, we now determine those 
combinations which form other bases. 



Triples of MU bases in C 4 

To begin, consider the MU vectors hi, ... , /14, in Eqs. (|15ll6p which exist for all values of 
x G [0, 7r]. Calculating their inner products, one finds that they only form an orthonormal 
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Figure 2: Plot of Eqs. (|13j) in the £+-plane implying that, for cos their solutions are 

given by = ±isina. 



basis of C 4 if y = y' and z = z' . Thus, for each value of x, the pair {J, F^x)} may be 
complemented by a third MU basis taken from the two-parameter family 



H4(y,z) 



( 1 
1 



iy 



1 
1 

e iy 
- P iy 



1 



1 \ 
-1 



(24) 



In other words, there is a three parameter-family of triplets of MU bases {J, F^(x), H^(y, z)} 
in dimension d = 4. This family was not known before. 

If x = 7r/2, additional MU vectors ji, ■ ■ ■ ,J4, and ki,...,ki, have been identified, cf. 
Eqs. (|21H23|) . Calculating the scalar products within each group, one sees that two further 
orthonormal two-parameter bases emerge, 



■h(r,s) 



K 4 (t,u) 
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e ir 
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(25) 



(26) 



9 



if the conditions r = r' , s = s', and t = t' , u = u', respectively, are satisfied. No other com- 
binations of the MU vectors can form inequivalent bases so that the matrices in Eqs. (|24H26|) 
represent all possible choices of a MU basis. Permuting appropriate rows and columns of 
the matrices J4 and K4 transforms them into H4; thus, the triples {/, F±(tt/2), J±(r, s)} 
and {I, F^(tt/2), K±(t, u)} are equivalent to {I, i*4(7r/2), Hi(y, z)}. 



Quadruples and quintuples of MU bases in C 4 

Let us begin by noting that sets of four MU bases cannot exist away from x 
two matrices H^(y,z) and H^y 1 ,z') are MU since 



74(2/^1 (*/) = h\(y)h 2 (y' 



tt/2. No 



(27) 



only hold if 



1 ± e i{y - y,) 



1 



(28) 



however, these equations have no solution for any values of y and y' . A similar argument 
shows that there are no values of z and z' such that the matrices H^(y,z) and H^{y,z') 
are MU. 

We now show that for x = tt/2 the bases H^(y, z), J±(r, s) and K^it, u) give rise to four 
and five MU bases if the free parameters are chosen appropriately. An argument similar to 
the one just presented shows that no two bases within either the family J±{r, s) or K^t, u) 
are MU. Thus, any quadruple of MU bases must contain bases from different families. 



The inner products 



h\(y)h( 



h\(y)j2(r) , 4(y)j'i( r ) and h Uv)h(r 



ulus 1/2 if there are values for y and r such that the equations 



1 + e ir ± (e~ iy + e i{r ~ y) 
l-e ir ± (e 



W _ e i{r-y) 



have mod- 

(29) 
(30) 



hold simultaneously. Upon introducing a factor of e ir / 2 ? Eqs. f)29[) are equivalent to the 
constraints 

' ' (31) 



r 

cos — 
2 



1 ± e~ % y\ 



r 

sin — 
2 



Consequently, one must have r = tt/2 , and thus e~ ty = ±i or y = tt/2 since < r, y < tt. 
An entirely analogous argument restricts the values of s and z: checking the inner products 
I ^1 ( r )^J3 (^) I ) ^lM e t c - t eus us that the matrices H 4 (y, z) and Ji(r, s) are mutually 

unbiased only if y = z = r = s = tt/2. We also find that the pairs { Ji(r, s), K^(t, u)} and 
{K^(t, u), H^x, y)} are MU only when all six parameters take the value tt/2. 
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We are now in the position to list all possible sets of MU bases in C 4 beyond {/, F^x)}, 

{I,F A (x),H A (y,z)}, 
{/, F 4 (tt/2), H^tt/2, tt/2), J 4 (tt/2, vr/2)} , (32) 
{/, F 4 (tt/2), F 4 (tt/2, vr/2), J 4 (vr/2, tt/2), K 4 {tt/2, tt/2)} . 

There is one three-parameter family of triples consisting of the one-parameter Fourier fam- 
ily F 4 (x) combined with two-parameter set H 4 (y, z); neither J 4 (r, s) nor K 4 (t, u) give rise to 
other triples since each of these sets of Hadamard matrices is equivalent to {/, F(n/2), H 4 (y, z)}. 
The three-dimensional set (|32p of MU bases in dimension d = 4 may be visualized as a 
cuboid defined by < x,y < ir and < z < tt/2. The reduction in the parameter 
range of z is due to the equivalence {/, F 4 (x), H 4 (y, z)} ~ {/, F 4 (tt — x), H 4 {tt — y,ir — z)} 
which follows from an overall complex conjugation. Each of the points in the cuboid corre- 
sponds to one triple while both the quadruple and the quintuple are located at the point, 
x = y = z = tt/2. 

Only one set of four MU bases exists, {I, F 4 (tt/2), i/ 4 (7r/2, n/2), J 4 (-7r/2, vr/2)}, since 
the other two candidates obtained by combining K 4 {ir/2, tt/2) with either J 4 (7r/2, tt/2) 
or i? 4 (7r/2, tt/2) are permutations of this quadruple. Finally, there is a unique way to a 
construct five MU bases which is easily seen to be equivalent to the standard construction 
of a complete set of MU bases in dimension four. 



4 Dimension d = 5 



As in dimensions two and three, there is a unique choice of a (5 x 5) dephased complex 
Hadamard matrix |18j . 
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7! 
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i i 
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UJ 
UJ^ 
UJ 2 



1 \ 

UJ 4 

UJ 2 
UJ ) 



(33) 



equal to the discrete (5 x 5) Fourier matrix, with uj = exp(27ri/5) denoting a fifth root of 
unity. The uniqueness of F§ is obtained by an analytical method that is far from elementary 

m- 

We have not found an elementary method to obtain a list of all vectors which are MU 
to the Fourier matrix F§. Instead, we will rely on earlier work [16] where those vectors 
have been constructed analytically by means of a computer program. 
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4.1 Constructing vectors MU to F 5 

The vector v = (l,e iai , . . . ,e ia4 )/y/5 G C 5 is MU to F 5 if it satisfies the conditions 



3=0 



y/E, k = 0..A, 



(34) 



defining ao = 0. According to [16], the solutions of these equations give rise to 20 vectors 
which can be arranged in four MU bases, 
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.(35) 



To obtain this result, Eqs. (134]) have been expressed as a set of coupled quadratic poly- 
nomials in eight real variables. Using an implementation [19] of Buchberger's algorithm 
[20} |2T] on the computer program Maple [22J , a Grobner basis of these equations has been 
constructed which leads to the 20 vectors given by the columns of the four Hadamard 
matrices above. It is important to note that no other solutions of Eqs. (I34p exist, a result 
which does not follow from the known methods to construct a complete set of six MU bases 
in C 5 . 

Each of the four matrices in (|35p is related to the Fourier matrix in a remarkably simple 
manner. In analogy to the unitary diagonal matrix used in Eq. (llOj) . define a diagonal 
unitary matrix 



D = diag(l, Co 1 , oj , ui 
(i) 



OJ 



with entries given by the first column of Hi and you find that 



(k) 



D k F, 



k = l,...,4. 



(36) 



(37) 



Using this observation, we can express the unique complete set of six MU bases for dimen- 
sion d = 5 as follows 



{/, F 5 ,HP,..., H^ 4) } = {/, F 5 ,DF 5 ,D 2 F 5 , D 3 F 5 , D 4 F 5 } , 
which will be useful later on. 



(38) 



12 



Next, we proceed to classify all smaller sets of MU bases of C 5 by combining subsets 
of the four Hadamard matrices He- in (|35p with the pair {/, F5}. For clarity, we now list 
the set of inequivalent classes which we will obtain. In addition to the pair {I, F5} and the 
complete set given in (|38[) there are two inequivalent triples as well as one quadruple and 
one quintuple: 

{I,F 5 ,H^},{I,F 5 ,H^}, 

{I,F 5 ,hP,H^}, (39) 

Triples of MU bases in C 5 

Select one of the four matrices given in ()35[) and adjoin it to the pair {/, i 7 ^}. You obtain 
four triples of MU bases with two immediate equivalences, namely, 

{I,F 5 ,H^} ee {I,F B ,DF S } ~ {I,F 5 ,D*F 5 } = {I,F B ,H^} (40) 

on the one hand, and 

{I,F 5 ,H^} ee {I,F 5 ,D 2 F 5 } ~ {I,F 5 ,D 3 F 5 } = {I,F 5 ,H®} (41) 

on the other. The equivalence ()40|) follows from multiplying the set {I, F B , DF§} with D 4 
from the left, rephasing the first basis with D from the right, using D 5 = / and swapping 
the last two matrices. A similar argument establishes the equivalence (|4ip . using D 3 instead 
of D 4 . 

Thus, it remains to check whether the triples T (1) = {I,F 5 ,H^} and T (2) ee {I,F 5 , 
(2) 

} are equivalent to each other. It turns out that these two triples are, in fact, inequiv- 
alent. More explicitly, this means that no unitary matrix U and no monomial matrices 
Mo, Mi and M2 can be found which would map into T^ 2 ' according to 

{/, F 5 , H^} -> {UIM , UF 5 M X , UH^M 2 } . (42) 

A proof of this statement is given in Appendix B. 

Quadruples of MU bases in C 5 

There are six possibilities to form quadruples by selecting two of the four matrices in Eq. 
(f35|) and adding them to the pair {I,F 5 }, Recalling that = D k F 5 , we identify the 
following equivalences which relate three quadruples each, 

{I,F 5 ,DF 5 ,D 2 F 5 } ~ {I,F 5 ,D 3 F 5 ,D*F 5 } ~ {/, F 5 , DF 5 , D*F 5 } , (43) 
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and 

{I,F 5 ,DF 5 ,D 3 F 5 } ~ {I,F 5 ,D 2 F 5 ,D 4 F 5 } ~ {/, F 5> D 2 F 5 , D 3 F 5 } . (44) 

To show the first equivalence in Eq. ([43 j) . for example, multiply its left-hand-side with D 3 
from the left, use the identity L> 5 = 1 and rearrange the bases appropriately. The other 
equivalences follow from analogous arguments. Thus, there are at most two inequivalent 
sets of four MU bases in C 5 , with representatives {/, F 5 ,H^,H^} and {I, F 5 ,H^,H^}, 
say. 

Interestingly, these two classes of MU bases are equivalent to each other leaving us with a 
single equivalence class of quadruples in dimension five, with representative {/, F$, }, 
say. To show this equivalence, we multiply the first quadruple with the adjoint of F§ from 
the left 

F B ,Hjp,H®} ~ {/, F 8 , FtffW, F^} . (45) 

using the identity F^ = FP, with some permutation matrix P, and swapping the first two 
bases. The action of F$ on the other two elements is surprisingly simple: the Hadamard 
matrix h9~ is mapped to itself, 



FiH, 



(46) 



(2) (3) 

up to a monomial matrix M, while iig is sent to H§ , 



FlH® =Hf>M>, 



(3) 



(47) 



again up to some monomial matrix M' . Both relations simply follow from working out the 
product on the left and factoring the result, e.g. 



(2) 



1 

71 



/ *(1) 8(2) 

w 3 s(l) u 2 s(2) 

lu 2 s(1) s(2) 

u 2 s(l) u A s(2) 

\ lu 3 s(1) lo a s(2) 



s(3) s(4) 8 (5) \ 

us(3) a (4) oj 4 s(5) 

w 3 s(3) los(4) lo a s(5) 

ws(3) w 3 s(4) s(5) 

s(3) wa(4) w 2 s(5) / 



fl -(3)_ D (2) P) (48) 



where the A;* 71 entry of the diagonal matrix is given by the sum of the k th column of 
denoted by s(k) = X/i-^ifc ' and P permutes the columns. Using these identities in 



Eq. (HSJ we find that 

{i,F s ,i#\HjP}~{i,F s M 1 \ni* i } 

the two quadruples are equivalent. 



(49) 
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Quintuples of MU bases in C 5 

Four sets of MU bases can be obtained by adding any three of the four matrices in Eq. 
(|35|) to the pair {J, F5}. It is not difficult to show that the four resulting sets of quintuples 
are equivalent to each other. Thus, there is effectively only one possibility to choose five 
MU bases in C 5 , with representative {I, F 5 , Hf \ Hf } }. 

Let us show now that this representative, which has been obtained by leaving out , 
is equivalent to the set {/, F5, H^}, for example. Indeed, the equivalence 

{/, F 5 ,DF 5 ,D 2 F 5 ,D 3 F 5 } ~ {/, F 5 , DF 5 , D 2 F 5 ,D 4 F 5 } , (50) 

follows immediately from multiplying the second set by D from the left and using D 5 = I, 

{I, F 5 ,DF 5 ,D 3 F 5 , D 4 F 5 } ~ {/, DF 5 , D 2 F 5 , D 3 F 5 , F 5 } . (51) 

Reordering the set of five matrices on the right reveals the desired equivalence with the 
quintuple {/, F5, flg , H^}. Effectively, the four matrices different from / undergo a 
cyclic shift under multiplication with D, and the remaining equivalences follow from shifts 
induced by D 2 and -D 3 , respectively. 

5 Summary and Discussion 

We have constructed all inequivalent sets of mutually unbiased bases in dimension two to 
five. Our approach is based on the fact that all complex Hadamard matrices are known in 
these dimensions. For dimensions up to d = 4, elementary arguments suffice to classify the 
existing sets of MU bases while dimension five requires some analytic results which have 
been found earlier using algebraic computer software. 



d 


2 


3 


4 


5 


6 


pairs 


1 


1 


oo 1 


1 


> oo 3 


triples 


1 


1 


oo 3 


2 


> oo 2 


quadruples 




1 


1 


1 


? 


quintuples 






1 


1 


? 


sextuples 








1 


? 



Table 1: The number of inequivalent MU bases for dimensions two to six where oo fc denotes 
a ^-parameter set; see text for details. 

The first four columns of Table Q] summarize the results obtained in this paper. All 
pairs of MU bases in dimensions two to five are listed in the first row, effectively reflecting 
the known classification of inequivalent Hadamard matrices; a continuous (one-parameter) 
set of inequivalent MU pairs only exists in dimension four. 
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The main results concern triples of MU bases in dimension four where we find a three- 
parameter family and in dimension five where we obtain two inequivalent triples. 

Finally, we have shown that there is only one class of both MU quadruples and MU 
quintuples in dimensions four and five. In all dimensions considered, there is a unique 
(i-tuple which can be extended to a complete set of (d + 1) MU bases using a construction 
presented in [23] , 

The last column of Table Q] contrasts these results with dimension six where the clas- 
sification of all complex Hadamard matrices is not known to be complete. The first entry 
shows that there is a three-parameter family of pairs of MU bases [23] (it has been con- 
jectured that the parameter space has, in fact, four dimensions |26j). Furthermore, it is 
possible to construct a two-parameter family of triples [25] using an idea taken from |27| . 
There is strong numerical |28[ [29] and analytical |16[ [T7] evidence to suggest that triples 
are the largest sets of MU bases in dimension six but this problem remains open. 

The notion of equivalence used in this paper (see Appendix \X§ is mathematical in 
nature; it captures all possible operations that leave invariant the conditions (PQ) for two 
bases to be mutually unbiased. Motivated by experiments, there is a finer equivalence of 
complete sets of MU bases based on the entanglement structure of the states contained 
in each basis |30[ 131] . For dimensions that are a power of two, a complete set of MU 
bases can be realized using Pauli operators acting on each two-dimensional subsystem. 
Two sets of MU bases are then called equivalent when they can be factored into the same 
number of subsystems. For d = 2, 4 this notion of equivalence also leads to a unique set 
of (d + 1) MU bases. However, for d = 8, 16, . . . complete sets of MU bases can have 
different entanglement structures even though they are equivalent up to an overall unitary 
transformation [301 EI] ■ 

The traditional approach to find complete sets of MU bases in prime-power dimensions 
via the Heisenberg-Weyl group or by using finite fields is constructive and, therefore, does 
not exclude the existence of other inequivalent complete sets. The present approach is, in 
contrast, exhaustive: we are able to affirm that the known complete sets for 2 < d < 5 are 
unique (up to equivalence). Their uniqueness has been shown earlier for d < 4 [32] while 
|33j contains a proof for 2 < d < 5 in a Lie algebraic setting. We find it appealing that 
it is possible to prove the uniqueness of complete sets of MU bases in low dimensions by 
elementary methods. 
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A Equivalent sets of MU bases 

Many sets of MU bases are identical to each other. To simplify the enumeration of all sets 
of MU bases we introduce equivalence classes and a standard form of sets of MU bases. 

Each set of (r + 1) MU bases in C d corresponds to a list of (r + 1) (with r < d) 
complex matrices H p , p = 0, 1, . . . , r of size (d x d). Two such lists {Hq, Hi, ... , H r } and 
{H'q, H[, . . . , H' r } are equivalent to each other, 

{Hq,Hi, . . . ,H r } ~ {H' ,H l , . . . ,H r } (52) 
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if they can be transformed into each other by a succession of the following four transfor- 
mations: 

1. an overall unitary transformation U applied from the left, 

{H , H u . . . , H r } -> U{H , H u ..., H r } = {UH , UH U UH r } , (53) 
which leaves invariant the value of all scalar products; 

2. (r+1) diagonal unitary transformations D p from the right which attach phase factors 
to each column of the (r + 1) matrices, 

{H , H 1 ,...,H r }^ {H D , H X D X , H r D r } ; (54) 

these transformations exploit the fact that the overall phase of a quantum state drops 
out of from the conditions of MU bases; 

3. (r + 1) permutations of the elements within each basis, 

{H , H 1 ,...,H r }^ {H P , ffiPi, • . . , H r P r } , (55) 

which amount to relabeling the elements within each basis by means of unitary per- 
mutation matrices P n satisfying PP T = I; 

4. pairwise exchanges of two bases, 

{• • • ,H p , . . . ,H p i, ...}->{.. .,H p i, . . .,H p , ...} , (56) 
which amounts to relabeling the bases. 

5. an overall complex conjugation 

{H , H 1 ,...,H r }-^{H^,Ht,..., H* r } (57) 
which leaves the values of all scalar products invariant. 

These equivalence relations allow us to dephase a given set of MU bases. The resulting 
standard form {I, Hi, ... , H r } is characterized by four properties: (i) the first basis is 
chosen to be the standard basis of C d described by Hq = I, where I is the (d x d) identity 
matrix; (ii) the remaining bases are described by (complex) Hadamard matrices: each of 
their matrix elements has modulus \j\fd; (iii) the components of the first column of the 
matrix H\ are given by 1/ \[d; (iv) the first row of each of the Hadamard matrices H\ to 
H r has entries 1 / \fd only. 
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B Inequivalent triples of MU bases in C 5 



We show that the two classes of triples of MU bases given by = {I,F B ,H^ } and 
7~( 2 ) = |J ) p 5j JJg } are inequivalent. In a first step, we explain that it is sufficient to 
search for equivalence transformations generated by matrices of a special form. In a second 
step we show that a contradiction arises if one assumes that the triples and are 
equivalent. 

Let us begin with a general remark about the structure of equivalence classes of sets 
of MU bases M. = {I, B%, . . . , B r } of C d for all r G {1, . . . , d — 1}. For convenience, we 
assume that the first basis equals the identity, i.e. the set is given in standard form. As 
explained in Appendix [A] all sets of MU bases equivalent to M are obtained as follows, 

M -»■ M' = {UM , UB\M\, UB T M r } ; (58) 

with a unitary U and (r + 1) monomial matrices Mj being a product of diagonal unitaries 
with permutation matrices; to keep the notation simple we do not reorder the (r + 1) bases 
within M' . For the set M! to be in standard form, one of the bases in Ai, say B p , must 
be mapped to the identity. As a consequence, the overall unitary transformation U must 
have a particular form, namely 

U = NBl, (59) 

where N is some monomial matrix and B p is one of the matrices contained in the set M. 
In view of Eq. (|59p we are lead to determine the action of Fl and (H^)^ on the triple 
as well as the action of and (H^)^ on the triple . It turns out that both 
triples are invariant under these global transformations as we have the equivalences 

NFlT {l) ~ NT (1) ~ N{H^ ] ^T W , (60) 

and 

AF 5 V (2) ~ AT (2) ~ A(tff)tT (2) . (61) 

The first equivalence in (|60p follows from using F^ = F^P and Eq. (|46p while the second 
one also requires the identity 

(H^F 5 = H^M, (62) 

with some monomial matrix M . The equivalences (161 p are derived in a similar way. 

Consequently, we can always remove the effect of the matrices Bp in the global trans- 
formations (|59p which leaves us with 

{/, F B ,H®} -> {A/M , NF 5 M h NH^M 2 } , j = l,2, (63) 
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where N, Mi and M 2 are monomial matrices, and up to rearranging terms. The non-zero 
entries of the monomial matrix iV must, in fact, be fifth roots of unity but we will not need 
this fact0. 

Using the restricted transformations shown in Eqs. (|63p . the triples {/, F5, Hg } and 

(2) 

{I,F^,H^ } are equivalent to each other only if either 

NF 5 = F5M1 and Hf ) M 2 = NH^ ] , (64) 

or 

NF 5 = F 5 (2) Mi and F 5 M 2 = NH { 5 1] , (65) 

hold for some monomial matrices Mi and M 2 . The choice Mo = iV -1 = N* in Eqs. f)63|) 
ensures that the identity will be mapped to the identity. 
Eqs. (|64p will now be shown to imply the identity 

AF 5 = F 5 M (66) 

for some monomial matrix M while A is a diagonal matrix with fifth roots of unity as 
nonzero entries, not proportional to the identity, A 7^ cl, c € C. However, Eq. (|66p only 
holds if A is a multiple of the identity. This contradiction implies that there are no matrices 
N,Mi,M2 such that Eqs. ([64"]) hold. Since Eqs. ([65]) also imply Eq. ([66]) with a (possibly 
different) diagonal matrix A 7^ cl, c E C, the triples and cannot be equivalent. 
Use i?5 = D^F^,j = 1, 2, to express the second equation in ([64"]) as 

D 2 F 5 M 2 = NDF 5 = NDN^NF 5 = DNF 5 , (67) 

introducing D = NDN^ = PDP T . Thus, the matrix D is obtained from D by reordering 
its diagonal elements according to the permutation P defined via N = PE, with some 
unitary diagonal matrix E. Combining this equation with the first one in (j64p leads to 
D 2 F 5 M 2 = DF 5 Mi, or 

D j D 2 F 5 = F 5 M X M\ (68) 

which is identical to (|66p upon defining A = D^D 2 and M = M\M\ which, as a product 
of two monomial matrices, is another monomial matrix. Since no permutation of the 
elements on the diagonal of D' = diag(l, w 4 , lo, uj, lo a ) produces the inverse of D 2 or a 
multiple thereof, we have A 7^ cl. Using the pair (j65l) instead of ([64]) also leads to an 
equation of the form (]66p with replaced by D which, however, cannot be a multiple of 
the inverse of D 2 , leading again to A / cl. 

1 Assume that N has a nonzero element different from a fifth root, say e la . This makes it impossible to 
transform T' 1 ' into standard form using right multiplication by monomial matrices unless the other nonzero 
elements of iV also equal e la . It follows that iV must be a permutation matrix P apart from a phase factor, 
N — e la P. Thus the matrices M p must have a common factor of e~ %a which, however, is irrelevant for the 
definition of MU bases 
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We now show that Eq. (|66p only holds if the matrix A is proportional to the identity. 
Write the monomial matrix M in (|66|) in the form 

M = PA" , (69) 

where P is a permutation matrix and A" is a diagonal matrix with entries having modulus 
one only. Denoting the inverse of A" by A', Eq. (|66p takes the form 

AF 5 A' = F 5 P. (70) 

Let us write A = diag(a,/3, . . . , e) with phase factors a,/3, etc, and similarly for A', and 
consider the simplest case P = I. Then the matrix relation ([70p reads explicitly 
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The conditions resulting from the first row immediately imply that the elements on the 
diagonal of A' are all equal to a*, or A' = a* I. The conditions of the first column imply 
that the matrix A is also a multiple of the identity, namely A = al. This contradicts the 
fact that the matrix A is different from a multiple of the identity. 

Let us now drop the restriction the P = I. The effect of P acting on F5 from the right 
is to permute its columns. The first row of -F5 will not change under this operation. Under 
the action of P, the first column will either stay where is is or it will be mapped to one 
of the four others. In the first case, we can immediately apply the argument given above 
to derive a contradiction. In the second case, it it straightforward to see that a similar 
argument still applies involving the first row of the matrices and that column which is the 
image of the first column. Thus, all possible choices of the monomial matrix M in (I66p 
require A to be a multiple of the identity — which it is not. 

Finally, we consider the action of an overall complex conjugation (I57p on either of the 
triples. We find that the set of three MU bases, remains invariant under complex 

conjugation 

r«* = {/,F 5 *,^ 1) *}~{/,F 5 ,^ 4) }~r«. (72) 

Similarly, complex conjugation maps to itself, T^* ~ T^. In summary, we have 
shown that the equivalence relations (|53p to f)57[) cannot transform the triple into T*- 2 -* 
or vice versa, i.e. these triples are inequivalent. 
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